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Abstract 

In this paper, we calculate properties of the spin polarized asymmetrical nuclear matter and 
neutron star matter, using the lowest order constrained variational (LOCV) method with the A Vis, 
Reid93, UV14 and AV14 potentials. According to our results, the spontaneous phase transition to 

m 

a ferromagnetic state in the asymmetrical nuclear matter as well as neutron star matter do not 
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I. INTRODUCTION 



The magnetic properties of polarized neutron matter, nuclear matter and neutron star 
matter have crucial role for studying the possible onset of ferromagnetic transition in the 
neutron star core. The most interesting and stimulating mechanisms that have been sug- 
gested for magnetic source of pulsars, is the possible existence of a phase transition to 
a ferromagnetic state at densities corresponding to the theoretically stable neutron stars 
and, therefore, of a ferromagnetic core in the liquid interior of such compact objects. Pul- 
sars are rapidly rotating neutron stars with strong surface magnetic fields in the range of 
10 12 - 10 13 Gauss Such a possibility has been studied by several authors using 

different theoretical approaches [4-29], but the results are still contradictory. In most calcu- 
lations, neutron star matter is approximated by pure neutron matter, as proposed just after 
the discovery of pulsars. Vidana et al. 21| and Vidana and Bombaci 22| have considered 
properties of spin polarized neutron matter and polarized asymmetrical nuclear matter us- 
ing the Brueckner-Hartree-Fock (BHF) approximation by employing three realistic nucleon- 



nucleon interactions, Nijmegen II, Reid93 and NSC97e respectively. Zuo et al. [25| have 
also obtained properties of spin polarized neutron and symmetric nuclear matter using same 
method with AVis potential. The results of those calculations show no indication of ferro- 
magnetic transition at any density for neutron and asymmetrical nuclear matter. Fantoni 
et al. 20j have calculated spin susceptibility of neutron matter using the Auxiliary Field 



Diffusion Monte Carlo (AFDMC) method employing the AU6 + UIX three-body potential, 
and have found that the magnetic susceptibility of neutron matter shows a strong reduction 
of about a factor 3 with respect to its Fermi gas value. Baldo et al. [26|], Akmal et al. 



27| and Engvik et al. 28( have considered properties of neutron matter with AVi$ po- 



tential using BHF approximation both for continuous choice (BHFC) and standard choice 
(BHFG), variational chain summation (VCS) method and lowest order Brueckner (LOB) 



v. On the other hand some calculations, like those of Brownell and Callaway 



4] 



5] considered a hard sphere gas model and showed that neutron matter becomes 



9J found that the inclusion 



respective 
and Rice 

ferromagnetic at k F ps 2.3/m~ 1 . Silverstein [8j and Ostgaard 
of long range attraction significantly increased the ferromagnetic transition density (e.g., 
Ostgaard predicted the transition to occur at kp ~ 4.1/m _1 using a simple central potential 
with hard core only for the singlet spin states). Clark [6] and Pearson and Saunier [h]] 



2 



calculated the magnetic susceptibility for low densities (kp < 2fm x ) using more realistic 



interactions. Pandharipande et al. 



Ill ], using the Reid soft-core potential, performed a vari- 



ational calculation arriving to the conclusion that such a transition was not to be expected 
for kp < 5/m _1 . Early calculations of the magnetic susceptibility within the Brueckner 
theory were performed by Backmann and Kallman 12) employing the Reid soft-core poten- 



tia 



and results from a correlated basis function calculation were obtained by Jackson et al. 
14] with the Reid v6 interaction. A different point of view was followed by Vidaurre et al. 
21] , who employed neutron- neutron effective interactions of the Skyrme type, finding the 
ferromagnetic transition at kp m 1.73 — 1.97/m . Marcos et al. 16( have also studied the 
spin stability of dense neutron matter within the relativistic Dirac-Hartree-Fock approxi- 
mation with an effective nucleon-meson Lagrangian, predicting the ferromagnetic transition 
at several times nuclear matter saturation density. This transition could have important 
consequences for the evolution of a protoneutron star, in particular for the spin correlations 
in the medium which do strongly affect the neutrino cross section and the neutrino mean 
free path inside the star [30]. 

In our pervious works, the properties of unpolarized asymmetrical nuclear matter have 
been considered by us u sing the lowest order constrained variational (LOCV) method 31] 
with the AV\s potential 32] . We have found that the energy per particle of nuclear matter 
alter monotonically and linearly by quadratic asymmetrical parameter. 



Recently, we have computed the properties of polarized nentron matter^ and polarized 
symmetrical nuclear matter [34[ such as total energy, magnetic susceptibility, pressure, etc 
using the microscopic calculations employing LOCV method with the AV\g potential. We 
have also concluded that the spontaneous phase transition to a ferromagnetic state in the 
neutron and asymmetrical nuclear matter does not occur. In this work, we intend to calculate 
the properties of spin polarized asymmetrical nuclear matter and neutron star matter using 
the LOCV method employing the AV 18 [32|, Reid93 [35|, UV U [Sjj and AV U [371 potentials. 



II. LOCV METHOD 



The LOCV method which was developed several years ago is a useful tool for the deter- 
mination of the properties of neutron, nuclear and asymmetric nuclear matter at zero and 
finite temperature. The LOCV method is a fully self-consistent formalism and it does not 
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bring any free parameters into calculation. It employs a normalization constraint to keep the 
higher order term as small as possible. The functional minimization procedure represents 
an enormous computational simplification over unconstrained methods that attempt to go 



beyond lowest order 



31 



38. 



39|]. 



In this method we consider a trial many-body wave function of the form 

4> = F<j>, (1) 

where <fi is the uncorrelated ground state wave function (simply the Slater determinant 
of plane waves) of A independent nucleon and F = F(l---A) is an appropriate A-body 
correlation operator which can be replaced by a Jastrow form i.e., 

F = SHf(ij), (2) 

i>j 

in which S is a symmetrizing operator. We consider a cluster expansion of the energy 
functional up to the two-body term, 

The one-body term E± is total kinetic energy of the system. The two-body energy E2 is 

E 2 = ^rE<V'Kl2)|v-J*>, (4) 



2^4 , 



where 



u(12) = -£[/(12),[V? 2 ,/(12)]] + /(12)V(12)/(12), /(12) and V(12) are the two-body 
correlation and potential. For the two-body correlation function, /(12), we consider the 



following form 



31 



3 

l 

k=i 

where, the operators 0^(12) are given by 



/(12) = E/ (fc) (12)0^(12), (5) 



(^-3) ( i2) = 1, ( 2 - + ls 12 ), (i-i-Su), (6) 

and 5*12 is the tensor operator. 

S12 = 3(cri.f)(cr 2 .f) - 0-1.0-2 (7) 
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Now, we can minimize the two-body energy Eq.(j4]), wit 
function fj- 1 ' but subject to the normalization constraint 



i respect to the variations in the 
3l|, 



TE^>l-f(12)Nj)a = 0, 

A a 



(8) 



where in the case of spin polarized neutron matter the function h$ z (r) is defined as 

-1/2 

/ l-.i /-•'. . \ \ " 

hsAr) 



1 9 /Jjfe) 



; S z = ±l 



l 



; S z = 



(9) 



here v is the degeneracy of the system. From the minimization of the two-body cluster 
energy, we get a set of coupled and uncoupled differential equations the same as presented 
in Ref. 311 ] . We can get correlation functions by solving the differential equations and in 
turn these functions lead to the two-body energy. 



III. SPIN POLARIZED ASYMMETRICAL NUCLEAR MATTER 



Spin polarized asymmetrical nuclear matter is an infinite system that is composed of spin 
up and spin down neutrons with densities pffl and pffl respectively, and spin up and spin 
down protons with densities pp and . The total densities for neutrons(p n ), protons 
(p p ), and nucleons (p) are given by: 



Pp = Pp ] + Pp\ Pn = P { n ] + Pn\ P = Pp + Pn 



(10) 



Labels 1 and 2 are used instead of spin up and spin down nucleons, respectively. One can 
use the following parameter to identify a given spin polarized state of the system, 



r Pp Pp tr Pn Pn 

Op , n 



(11) 



P P 

S p and S n are proton and neutron spin asymmetry parameters, respectively. Total polariza- 
tion defined as, 



5 = 5 n + 5 P 



(12) 



For the unpolarized case, we have 5 n = 5 P = 0. 
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Asymmetry parameter describes isospin asymmetry of the system and is denned as, 

p= Pn^Pp_ (13) 
P 

Pure neutron matter is totally asymmetrical nuclear matter with (3 = 1 and symmetrical 
nuclear matter has (3 = 0. In the case of spin polarized asymmetrical nuclear matter, the 
energy per particle given by 

E(p,5,(3) = E 1 (p,5,(3) + E 2 (p,6,py (14) 

where Ei, the kinetic energy contribution given by, 



MM/?) = l^f(l + / 3 + 25 n ) l + (i + /3-25 n 



+(l-/? + 2<y§ + (l-/?-2<y§}, (15) 

where kp = (3/2n 2 p) 1 ^ 3 is Fermi momentum of unpolarized symmetrical nuclear matter. In 
expression $Hj), the two body energy can be calculated by the semi-empirical mass formula 
3lL 



E 2 {p,5,f3) = (3 E 2neum + (1 - (3 )E 2snucm , (16) 

where the powers higher than quadratic are neglected. E 2neum is polarized pure neutron mat- 
ter potential energy, and E 2snucm is polarized symmetrical nuclear matter potential energy 
that both of them have been determined using the LOCV method employing a microscopic 



point of view 



33 



34]. 



The energy per particle of the polarized asymmetrical nuclear matter versus density at 
different values of the spin polarization and isospin asymmetry parameter have been shown 
in Figs. 1-4 for the AVi%, Reid93, UVu and AVu potentials, respectively. As it can be seen 
from these figures, the energy of polarized asymmetrical nuclear matter for various values 
of isospin asymmetry parameters become more repulsive by increasing the polarization for 
all relevant densities. There is no crossing of the energy curves of different polarizations. 
Reversibly, by increasing density, the difference between the energy of nuclear matter at 
different polarization becomes more sizable. Indeed, the ground state of the system is that 
of unpolarized matter in all ranges of density and isospin asymmetry considered. This 
behavior is common for all potentials used in our calculations. 
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Another quantities that help us to understand ferromagnetic phase transition, are mag- 
netic susceptibility and Landau parameter. In what follows, we want to consider these two 
quantities. The magnetic susceptibility, x> which characterizes the response of a system to 
the magnetic field and gives a measure of the energy required to produce a net spin alignment 
in the direction of the magnetic field, is defined by 

( dM \ 

* = UL' (17) 

where M is the magnetization of the system per unit volume and H is the magnetic field. 
We have calculated the magnetic susceptibility of the polarized asymmetrical nuclear matter 
in the ratio x/Xf form. By using the Eq. (1171) and some simplification, the ratio of \ to the 
magnetic susceptibility for a degenerate free Fermi gas xf can be written as 

X 2 E F 



(18) 



'5=0 



where Ep = h k F /2m is the Fermi energy. After doing some algebra the ratio x/Xf takes 
the form, 

-i 

(19) 



X_ 

Xf 



2 {j3 2 {2z(XF/x)neum ~ (XF / X) snucm] + (Xf/x) snucm + ( 1 ~ 2?) - l) 
1+ (l + /3) 5 /3 + (1-/3)5/3 



where (x/XF)neum and ixl Xf) amuan is magnetic susceptibility of pure neutron matter [331 ] . 
and asymmetry nuclear matter [3J], respectively. In Fig. 5, we have plotted the ratio x/Xf 
as a function of density at several values of the isospin asymmetry for the AVis, Reid93, 
UV14 and A V14 potentials. For all used potentials, we see that the value of the ratio x/Xf 
decreases monotonically with the density, even at high densities. It shows, there is no 
magnetic instability for any values of asymmetry parameter. If such an instability existed, 
the value of ratio x/Xf changes suddenly unlike previous treatment. 

As it is known, the Landau parameter, Go, describes the spin density fluctuation in the 
effective interaction. G is simply related to the magnetic susceptibility by the relation 

f F - TTg~ ' (20) 

where m* is the effective mass. A magnetic instability would require Go < — 1. Our results 
for the Landau parameter have been presented in the Fig. 6 for the AV^s, Reid93, UV\4 
and AVu potentials. It is seen that the value of G is always positive and monotonically 



increases by increasing the density. This shows that for all used potentials, the spontaneous 
phase transition to a ferromagnetic state in the asymmetrical nuclear matter does not occur. 

The equation of state of polarized asymmetrical nuclear matter, P(p,(3,5), can be ob- 
tained using 

P{p , M=p i a -E?im. (21) 

In Fig. 7, we have shown the pressure of asymmetrical nuclear matter as a function of 
density (p) at different polarizations for various choice of asymmetry parameter with 
the AVis potential. This figure shows that the equation of state becomes stiffer by increasing 
the polarization for all isospin asymmetry. 

IV. POLARIZED NEUTRON STAR MATTER 

Indeed the neutron star matter is charge neutral infinite system that is mixture of asym- 
metrical nuclear matter and leptons, specially electrons and muons. The energy per particle 
of polarized neutron star matter, E nsm can be written as, 

E nS m = E + Ei . (22) 

where E is the nucleonic energy contribution which given by Eq. (TH|) and Ei is leptonic 
energy contribution obtained as follows, 

Ei = E E M 2 + ^ 2 f /2 . (23) 

i=e, /i k<ki F ) 

After some simplification the above equation leads to 



1 ~ 8 . ^ (F) 3 
° i=e, ix r X\ ' 

>0 _ %uW 



xf\2xf )2 + l)Jxf )2 + 1 - sinh" 1 x\ F) 



(24) 



where x t = K h Mc . The conditions of charge neutrality and beta equilibrium impose 
the following constraints on the calculation of energy of neutron star matter 

Pn = Pp + He He = /V • (25) 

Pp = Pe + Pfi ■ (26) 
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In Fig. 8, we have shown the energy per particle of the neutron star matter at various values 
of spin polarization as a function of density for the AVis, Reid93, UVu and AVu potentials. 
The same as Figs. 1-4, this figure have also shown no phase transition to a ferromagnetic 
state. 

The treatment of magnetic susceptibility and landau parameter of neutron star matter 
versus density for different polarization, have been shown in Figs. 9 and 10, respectively. 
This treatments are the same as the treatments of the polarized neutron matter, symmet- 
rical and asymmetrical nuclear matter, and it does not show any magnetic instability for 
the neutron star matter. This behavior has been observed for all potentials used in our 
calculations. 

For the AVi$ potential, the pressure of neutron star matter have been presented as a 
function of density p at different polarizations in Fig. 11. We see that the equation of state 
of neutron star matter becomes stiffer by increasing the polarization. 

V. SUMMARY AND CONCLUSIONS 

The purpose of this paper was to calculate the properties of the spin polarized asym- 
metrical nuclear matter and neutron star matter employing the lowest order constrained 
variational technique with the AVi 8 , Reid93, UV U and AVu nucleon-nucleon potentials. 
After introducing the LOCV method briefly, we calculated the energy per particle of the 
polarized asymmetrical nuclear matter and showed that the force becomes more repulsive 
by increasing the polarization for all relevant densities. No crossing of the energy curves 
was observed and the ground state of the system was found to be that of the unpolarized 
matter in all ranges of densities and isospin asymmetry. The magnetic susceptibility was 
computed and shown to remain almost constant by increasing the asymmetry parameter. 
For the polarized neutron star matter, we showed that there is no magnetic instability and 
the equation of state becomes stiffer by increasing the polarization. 
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FIG. 1: The energy per particle of the polarized asymmetrical nuclear matter versus density (p) for 
different values of the spin polarization (5) with the AVis potential at (3 = 0.0 (a), j3 = 0.3 (b), 
(3 = 0.6 (c), 0= 1.0 (d). 
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FIG. 2: The energy per particle of the polarized asymmetrical nuclear matter versus density(p) for 
different values of the spin polarization (5) with the Reid93 potential at (3 = 0.0 (a), f3 = 0.3 (b), 
(3 = 0.6 (c), 0= 1.0 (d). 
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FIG. 3: The energy per particle of the polarized asymmetrical nuclear matter versus density(p) for 
different values of the spin polarization (5) with the UV14 potential at (3 = 0.0 (a), (3 = 0.3 (b), 
(3 = 0.6 (c), (3 = 1.0 (d). 
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FIG. 4: The energy per particle of the polarized asymmetrical nuclear matter versus density (p) for 
different values of the spin polarization (<5) with the AV14 potential at (3 = 0.0 (a), (3 = 0.3 (b), 
/3 = 0.6 (c), 0= 1.0 (d). 
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FIG. 5: The magnetic susceptibility of the polarized asymmetrical nuclear matter as the function 
of density (p) for different values of asymmetry (/?) with the AVis (a), Reid93 (b), UV14 (c) and 
AVu (d) potentials. 
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FIG. 6: The Landau parameter of the polarized asymmetrical nuclear matter as the function of 
density (p) for different values of asymmetry (/3) with the AVi$ (a) , Reid93 (b) , UVu (c) and AVu 
(d) potentials. 
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FIG. 7: The equation of state of polarized asymmetrical nuclear matter for different values of the 
spin polarization (<5) with the AV\% potential at j3 = 0.0 (a), (3 = 0.3 (b), = 0.6 (c), (3 = 1.0 (d). 



18 




0-1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 r 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

p(fm- 3 ) 




i -\ 1 1 ■ 1 ■ 1 1 1 ■ 1 ■ 1 1 1 1- 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

p(fm - 3 ) 



FIG. 8: The energy per particle of the polarized neutron star matter versus density (p) for different 
values of the spin polarization (S) with the A V±s (a), Reid93 (b), UV\± (c) and AV14 (d) potentials. 
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FIG. 9: The magnetic susceptibility of the polarized neutron star matter versus density (p) with 
the -AVi8, Reid93, UVu and AVu potentials. 




FIG. 10: The Landau parameter, Go, of polarized neutron star matter as function of density(p) 
with the AVi8, Reid93, UVu and AV14 potentials. 
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FIG. 11: The equation of state of polarized neutron star matter for different values of the spin 
polarization (<5) with the AVi% potential . 
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